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The frequency-dependent conductivity is studied for the one-dimensional Hub- 
bard model, using a selection rule, the Bethe ansatz, and symmetries associated with 
conservation laws. For densities where the system is metallic the absorption spec- 
trum has two contributions, a Drude peak at u = separated by a pseudogap from 
a broad absorption band whose lower edge is characterized by a non-classical critical 
exponent. Our findings shed new light on the "far infrared puzzle" and other optical 
properties of metallic organic chain compounds. 
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Organic chain materials [|TJ have a broken symmetry ground state (GS) and at low critical 
temperatures T tr undergo transitions to unusual T > T tr states. Recent T > T tr optical 
measurements over a huge frequency range |l|fj] established a pseudogap feature together 
with a zero- frequency mode for the (TMTSF) 2 X salts (TMTSF for tetramethyl tetraselena 
fulvalene.) These experiments have raised several important yet unresolved issues. The first 
U was related to the relevance of calculations based on the doped ID Hubbard model. 
Although these materials have finite interchain hopping integrals i& and t c with t a » tj, » 
t c 0, hopping becomes less effective as the temperature or the frequency is increased and 
the along-chain hopping t a dominates. 

In this Letter we study the optical conductivity of the ID Hubbard model, <J\(uj) = 
2tiD5{uj) + a[ e9 {uj), where a[ e9 (u)) is defined in Eq. (P below. Schulz used the Bethe- 
ansatz (BA) solution to calculate D and to derive the total intensity of finite-frequency 
transitions. Our study refers to all parameter space where that intensity is significant. 
Analytical calculations of a[ e9 (uj) have so far been restricted to the insulator half filling 
phase in the limit of large onsite repulsion U ||. In this case the optical gap occurs at 
Eopt = U — 4t a , followed by an absorption band extending up to U + 4t a Our study 

confirms the expectation |1J that calculations based on the doped ID Hubbard model are 
relevant for the theoretical understanding of recent optical experiments |],|2[] . In the metallic 
phase we find that the finite-frequency transitions are limited to a well-defined band above 
a pseudogap (which is actually smallest at half filling). The occurrence of such pseudogap 
agrees with predictions from exact diagonalizations which have found practically almost 
no low-frequency (u 7^ 0) absorption at the metallic phase 0. By combining symmetries 
associated with conservation laws || with the approach of Ref . , we are able to derive the 
exponent for the frequency dependence above the absorption edge (and below the top edge 
for half filling). We confirm [|I| that the existence of a pseudogap in the charge excitations 
with the absence of a gap for spin excitations is an indication of charge-spin separation in 
the metallic state and that that separation is distinct from that of a ID Tomanaga-Luttinger 
liquid (TLL), in which both excitations are gapless 0. We characterize the nature of the 
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unusual metallic state which following Ref. [|TJ we call doped- Mott- Hubbard liquid (DMHL). 
Our study reveals that in the case of the ID Hubbard model previous naive TLL theoretical 
studies |L0| of a[ e9 (co) do not apply to values of ui larger than the pseudogap. In contrast to 



previous predictions || , we find that the phenomenon of suppressed far infrared conductivity 
is of ID character and follows from a parity selection rule which also occurs at higher 
dimensions. However, only at ID is this selection rule effective in preventing far infrared 
conductivity transitions. For instance, in the 2D case doping away from half filling induces 
significant mid-infrared absorption |7,|Tl|. These results seem to indicate that although for 



low temperatures just above T tr the small interchain coupling might lead to 2D character 
for some of the properties of the (TMTSF) 2 X salts, in what the conductivity along the 
chains is concerned the main role of that coupling is equivalent to small doping on the single 
chains, which otherwise remain of ID character, with the electrons confined in each chain. 

The ID Hubbard model in a chemical potential fi describes N interacting electrons 
and can be written as H = H SO (a) ~ n[N a — N], where Hso{4) = T + U[D — N/2], T = 
—taJ2j,a[ c ]<T C j+iv + h.c] is the "kinetic energy", D = Sj%,T%,J. measures the number of 
doubly occupied sites, N = 5Zj> C ] CT an d c Jcr are electron operators of spin projection 
a at site j = 1, ...,N a , and h^ a = c^ a Cj a . We choose a density n = N/N a in the interval 
< n < 1 with even N and zero magnetization and define hp = im/2. We use units such 
that — e = h — 1, where — e is the electronic charge. The conductivity can be written as 



TT _ |(HJ|0)| 2 



(J T g { U} ) = TrY. 5(w-w Vt o), (1) 



and is also given by a[ e9 (u) oc linu^o ulrnx p^ k ^ ; where J = —it a J2j,a 



is the current operator, the summation runs over energy eigenstates, uj v $ = E u — E is the 
excitation energy, xp{k, u) = —J2u^o i s the charge-charge response function, 

and n(k) = J2k',a c l'+k a c k',a- The formula (|1|) applies if the GS |0) is non-degenerate. This 
is true for an even number of sites N a if periodic (anti-periodic) boundary conditions are 
used for odd (even) values of y, respectively. The GS is then necessarily an eigenstate of the 
parity operator, P n , which moves electrons from sites j to N a + 1 — j, j = 1, N a , and has 



eigenvalues ±1. Importantly, in the present model P n commutes with the Hamiltonian and 
anticommutes with the current operator. This implies immediately that the states |0) and 
J|0) have opposite parities, a key selection rule for optical transitions. The final states \i>) can 



be characterized in terms of holons, antiholons, spinons jl^JTql , and a charge-transfer band 
(CTB) |14j]. The holon/antiholon and spinon bands describe low-energy charge and spin 
excitations, whereas the CTB is associated with the upper Hubbard band 0. We use the 
labels a = c, s,t for the holon/antiholon band, the spinon band, and the CTB, respectively 
n|, and the quantum number (3 for distinguishing between holons (j3 = — |) and antiholons 



(P = +|) |T3| . |T4] ] . In the present context, the excitations that couple to the GS are such that 

the s band is empty of spinons, the c band can be populated by holons and by zero or one 

antiholon, and the t band can have occupancy zero or one. The momentum variable has the 

form qj = j^-Ij 1 , where Jj* are successive integers or half-odd integers. In contrast to the case 

of electronic bands, the number of available momenta N* can change with band occupancies. 

These numbers are N* = N a , N* = N/2 - N t , and = N a — N + N t , where N t is the 

number of occupied momenta in the CTB, and the momentum band widths are Aq c = 2tt, 

Aq s = 2k F - 2irN t /N a , and Aq t = 2vr - Ak F + 2irN t /N a . The numbers N h _ L (holons), 

N h i (antiholons), and A^ t are good quantum numbers which obey the sum rules ||14|| , 
c,+ 2 — 

N a — N = —2 E/9=±i P N c,/3 = N c - 2N t, where N% = T.p = ±i N^. The GS is characterized 
by N h i = N t = and N h _i = N a — N, with a symmetrical holon occupancy in the c band 

C,+ 2 C, 2 

for momenta 2kp < \q\ < 7T. The GS a pseudo-Fermi momenta read q Fa = KN a /N a where 
N a = — N% and thus q Fc = 2k F , q Fs = k F , and q Ft = 0. The energy bands e Q a (q) are 
defined by Eqs. (110)-(112) and (C7)-(C9) of Ref. 0. [See also Figs. 7 and 8 of Ref. gH 
for e a (q) = e Q a (q) — e°(?i? a ) an d ck = c, s] The velocities v a (q) = de° a (q)/dq and v a = v a (q Fa ) 
and the mass 1/m* = \dv a (q) / dq\ q = qFa are also important quantities. In the parameter 
region of appreciable oscillator strength for optical transitions, which corresponds to n — 1 
and to densities close to n = 1 for large enough (but not too large) values of U 0] , the CTB 
width Wt = e? (0) — e°(7r — 2k F ) = e°(0) either vanishes (n = 1) or is very small. There is 
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nearly no oscillator strength for U < U*, where U* = U*(n) is such that E c (n, U*) = 0. Here 
E c (n, U) = E C = -2e°(vr) and E c > (and E c < 0) for U > U* (U < U*) and E c -+ U - it 



(E c — > E M h) as n — > (n — > 1). (E MH is the n — 1 Mott-Hubbard gap [[14]].) U* is a 
decreasing function of n given by [/* = (U* = it) for n — > 1 (n — > 0). 

The main transitions contributing to a\ e9 {uo) are of two types: (a) those leaving the 
band filling N h _ 11 N i, and iV t unchanged and (b) those changing these numbers by 

C ' 2 c '+2 

AAf fe _i = AAf ft i = AN t = 1. The former start, in principle, at = 0, while the latter 

c, 2 c,+ 2 

have an onset at E opt = W t — 2e®(2kp) > 0. We find below that for n ^ 1 (or n = 1) 
E opt is a pseudogap (or a gap). Importantly, for uj < E opt the final-state Hilbert subspace 
is spanned by states with no occupancy in the t band. In this case the excitations are 
described by holons and spinons and at low values of uj and n ^ 1 the quantum liquid is 
a TLL. On the other hand, for uj > E opt the final-state subspace is spanned by states also 
involving t occupation, the charge and spin excitations remaining separated and described 
by holons and spinons. For 1 this c, s, and t quantum liquid describes the DMHL. The 
collisions involving holons, antiholons, spinons, and t particles do not lead to energy and 
momentum transfer and only give rise to shifts in the phases of these elementary particles. 
We find below that for both the TLL and DMHL the critical exponent expressions are simply 
a linear superposition of the phase-shift parameters £ J aa , = 5 a ^ + Y^i=±i^^aa'(QFa,lQFa') 
with j = 0, 1. Here <& aa '{q, Q 1 ) are the two-particle forward-scattering phase shifts defined by 
Eqs. (103) and (B30)-(B38) of Ref. [[14]]. For holons and spinons the 8 parameters £ 3 aa , read 



e cc = 2&S = l/e = V^p, & = -C C = V& = VV2, and & = e cs = where K p is the 
TLL parameter defined in Ref. and ^\ a , are the entries of the transverse of the conformal- 



field theory (CFT) dressed-charge matrix [16|. The TLL critical expressions involve v c , v s , 



and K p . On the other hand, the low-energy (u — E opt ) > DMHL critical theory (above 
the pseudogap) is a three-component quantum liquid of the general type studied in Ref. 
and involves v c , v s , mjf, and the 18 parameters ^ aa i- Importantly, some of these parameters 
are n and U dependent and cannot be expressed in terms of the TLL quantity K p . (The 
n = 1 and (u — Emh) > critical theory is also of that type and involves v s , m*, and ml.) 



The transitions of type (a) do not exist at n = 1 and we show first that, as a consequence 
of the selection rule presented above, they have a very small weight for n < 1. These are 
excitations within the holon band which can be characterized in terms of single and multiple 
electron-hole excitations |T7|]. Single electron-hole excitations give no contributions to the 



finite-frequency absorption, while multiple excitations are expected to decrease very rapidly 
in intensity. For low to we evaluate a[ e9 (u) from Xp(k,to) ex / dx J dte^ kx ~^x P {x,t), with 
the charge-charge correlation function x P (x, t) obtained directly from CFT ||16|| . Double 
excitations, which would give a contribution a\(uS) = C3U 3 , are forbidden by the parity 
selection rule and thus C3 = 0, since these transitions require symmetrical changes of holon 
occupancies at q and —q, leaving the parity unchanged. This explains why the a; 3 absorption 
was not observed in optical experiments 0. Higher-order low-energy processes would give 
rise to ai(u) = Cju 7 for n ^ 1 and to (J\{oS) = Aj uo^' 2Kp ~ l for n = 1/j with j = 2,3,4... 
(due to Umklapp processes). [These two exact low-a; CFT expressions can be obtained by 
multiplying the perturbative expressions u 3 and u^ 2Kp ~ 5 , respectively, of Ref. |L(J by the 
factor cj 4 , which is missing in the latter two expressions.] Since these processes require mul- 
tiple excitations in the holon band, they have extremely low spectral weight and C7 , Aj ~ 0. 
These considerations explain why numerical studies yield so weak features at low frequen- 
cies [0. While in ID the P n eigenvalues ±1 are directly related to the right and left state 
occupancies and the selection rule is very restrictive, at higher dimensions it is less effec- 
tive because there are many angular-momentum channels which lead to significant spectral 



weight at small values of 00 |TT|. The weakness of type (a) transitions implies that the onset 
of type (b) transitions at E opt represents a pseudogap which for n < 1 and U > U* looks like 
a true gap in an actual experiment [0,0. E opt , shown in Fig. 1 (a), increases with increasing 
U and decreases with increasing density. For n — 1 it coincides with the Mott-Hubbard gap, 
Emh — E C (1,U). Since larger interchain transfer integral tj, means larger effective doping 
for the single chains, this theoretical result is fully consistent with the experimental data of 
Fig. 3 of Ref. |IJ . The important processes for the type (b) transitions are those where in 
addition to a particle at q' in the t band a holon-antiholon pair is created in the c band at 
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q and q' — q (with q' = for n = 1). This leads to final states whose excitation energy has 
its minimum and maximum values at u — E opt and uj = E opt + Auj opt , respectively. Here 
Aoj opt = 8t a — 2ep where e h F = e°(7r) — e°(2/c^) = at n = 1. This spectral range agrees with 
large U expansions ||] and numerical calculations [^,0 . For n — 1 there are only these final 
states and for n < 1 the missing multipair processes contain nearly no spectral weight. The 
current operator cannot produce more than one doubly occupied site. Since from the results 



of Refs. ]14] ) |17|] we find that creating one t particle involves processes which create one of 
these sites, there will be almost no spectral weight above the upper edge at uj = E opt + Au opt . 

The evaluation of the line shape is a very difficult problem. We could solve it both for 
n < 1 and the region immediately above threshold and for n = 1 and energies immediately 
below the top absorption edge. Since the number operators N a are conservation laws || 
and the t band is quadratic for small q, the low-energy spectrum of the Conservation-Laws 
Hamiltonian, Hcl = Hso{&) — Ha eQ a (.lFa)N a — fiN a , is of the type studied in Ref. 0. The 
use of the non-linear band approach of that reference allows the evaluation of asymptotic 
expressions for the charge-charge correlation function of Hcl, Xp L ( x jt)- However, that 
approach does not provide finite-energy correlation-function expressions. Fortunately, the 
following symmetries and properties of H and Hcl allow solution of this problem [Q: (i) 
the leading term in Xp \ x >t) nas contributions only from transitions to a Hilbert subspace 
with fixed N a numbers; (ii) [H, Hcl] — 0, and thus H and Hcl have the same eigenstates; 
and (iii) for the final subspace relevant to the conductivity, the difference in energies is 
E op t = {i>\H\il)) — (iP\Hcl\?P) and is the same for all final states \ip) in that subspace. 
Using these results, we can show that the corresponding leading term in the Hubbard model 
correlation function is of the form xpi x ,t) = e lEoptt x^ L {x,t). This provides Xp{^i u ) = 
J dx J dte^ kx ~( u ~ Eovt ^x^ L ( x ,t) fc> r small values of both k and {uj — E opt ) > which for 
n < 1 leads to a[ eg (uj) = C( ^f° p *0 at the onset and a critical exponent £ = — | + 
|(\ZtI~ ~ £ct) +l(^st) • A s expected, this DMHL exponent cannot be expressed in terms 
of the TLL parameter K p only. It is shown in Fig. 1 (b) and approaches | both as n — ► 1 
at finite U and as U — > oo, and it increases with decreasing values of U. C increases with 



increasing n and for U < U* is vanishing small. Consistently with the n < 1 results, at 

n = 1 we find a r i 9 {uj) = C ( ) ^ for all finite values of U, in contrast to the naive 
\— 1/2 

— Emh) prediction of Ref. At n = 1 we could also find that a[ eg (cu) has a 

maximum at u; = u Q with £ Mfl < c^o < -Em// + 4t a and co> - > E MH + 4i a as U — > oo, and 
that for energies immediately below the top conductivity edge at -Em// + 8t a , <j[ e9 (o;) reads 
a[ e5, (u;) = (7 / ^ [ gAf g+ 8 ^ a ~"] ^ . At U large the ratio C'/C equals one and the conductivity 
absorption is symmetric around u = ujq = Emh + 4t a ~ U, in agreement with the results of 
Refs. ||,||. However, C'/C decreases with decreasing U and u is shifted to lower values of uj. 
In this case a\{uS) is not symmetric around ojq and a dependence <J\(u) oc u;~ 7 is expected for 
a conductivity descending region of intermediate uj values of the domain u G (u , E opt + 8t a ). 
(In contrast to the onset exponent (, 7 cannot be calculated within our critical theory and 
the use of TLL |@,10| to estimate it is an open question.) 



The X = C10 4 and T = 10K curves of Figs. 1 (B) of Ref. and Figs. 3 and 4 of Ref. @ 
correspond to a uj = Drude peak with a very small carrier density separated by a pseudogap 
E op t ~ 0.014 eV from an absorption band with 7 rs 1.3 and a width Auj opt ~ 1.00 eV. 
These -E OJ3< and Au opt values can be understood in the framework of the ID Hubbard model 
with effective values t a ~ 0.125eV [tight-binding model (TBM) calculations lead to 
t a ~ 0.250 eV because the energy width of the e c (q) band GS occupancy, 4t a — €p, is about 
twice that of the TBM band] and U/t a ~ 1.5 if the system has a nearly half-filled band, 
n ~ 0.995, in which case the optical threshold E opt « £^mh ~ 0.11 t a « 0.014 eV and Au; op t ~ 
8t a ~ 1.00 eV. Moreover, combining the theoretical values for Atu^ rs 8t a and the sumrule 
J^° 00 o'i(uj) 0, we find C ~ 1.2 <r[ eg (u> ) ~ 3600fi _1 cm _1 . A crucial text for our theoretical 
description is whether the exponent ( obtained from the observed sharp onset of absorption 
]T|fJ] agrees with the t a , U, and n values extracted from E opt and Auj opt . Except for the small 
onset temperature tail, we find from the above figures the linear behavior, In ( °k<, ^ ) ~ 
C ln( [<J ~ g ° pt] )+ln( r ) for fe^^I < 0.5, where C ~ 0.65 and C « 3530 fi^crrT 1 . The 
corresponding asymptotic theoretical conductivity curve is plotted in Fig. 2. Importantly, 
Fig. 1 (b) reveals that this exponent is consistent with the above U/t a and n values. 
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FIGURES 

FIG. 1. (a) The pseudogap E opt in units of t a and (b) the exponent C as function of n for 
different values of U. The solid (dashed) lines correspond to U > U* (U < U*). 

FIG. 2. The conductivity a[ e& '(u) of (TMTSF) 2 C10 4: for T = WK as a function of 
[uj— Eopt] /Eopt > [from Fig. 4 of Ref. Q] and the theoretical asymptotic curve c([ui— E op t]/ E pty" 
(solid line) with ( = 0.65 and C = 3530 Qr x cmr x . 
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